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DISTRIBUTION AND WAVE IRAQ FOR A LIMITED FAMILY 
OF TAPERED SWEPTBACK WINGS WITH SYMMETRICAL 
PARABOLIC -ARC SECTIONS AT ZERO LIFT 
By Julian H. Kainer 

SUMMARY 


The study of the supersonic pressure distribution and wave drag of 
sweptback wings at zero lift has been extended by means of linearized 
theory to include the combined effects of plan -form taper and curvature 
of profile. Calculations have been made for a specific wing with a 
supersonic leading edge and these results are applicable to a limited 
family of plan forms. Comparisons have been made of the pressure dis- 
tributions and the wave drag calculated by this exact linear solution 
with the results of an arbitrary simplified approximation based on 
formulas for untapered plan forms with symmetrical parabolic -arc sections. 
In addition, a comparison of the wing wave -drag coefficient has been made - 
with the results obtained from formulas for geometrically similar swept- 
back tapered wings with rhombic sections modified by a profile -correction 
factor . 


INTRODUCTION 


During the past few years, extensive investigations utilizing 
linearized theory have been in progress to predict the aerodynamic 
characteristics of a wide variety of wing configurations at supersonic 
speeds. The method of reference 1 was applied for the case of zero lift 
in references 2 and 3 to predict the wave-drag characteristics of 
untapered sweptback wings with symmetrical parabolic-arc sections nnfl 
in references 4 and 5 to predict the wave drag of tapered sweptback wings 
with rhombic profiles (symmetrical double-wedge sections with maximum 
thickness at the midchord) . Reference 6 considered the effect of 
chordwise location of maximum thickness on the supersonic wave drag of 
sweptback wings with symmetrical double -wedge profiles. Since many 
practical configurations Include combined plan-form taper and curvature 
of profile, it is of current interest to investigate this case. This 
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paper supplements the investigations of references 1 to 5 by presenting 
the calculations for the pressure distribution for a limited family of 
tapered sweptback plan forms with symmetrical parabolic-arc sections and 
supersonic leading edges at zero lift. 


SYMBOLS 


x, y, z Cartesian coordinates with origin located at wing apex 

tj coordinates of origin of source line 

V free-stream velocity in x-direction 

M free-stream Mach number 

0 X first partial derivative of disturbance-velocity 

potential in x-direction 


0xx> 0zz second partial derivatives of disturbance -velocity 

^ potential in x-, y-, and z-directions, respectively 

3 - \Jm 2 - 1 

A sweep of leading edge 


b 

c r 

ct 

X 


xi, yi 


c 


wing span 
root chord 


tip chord 

taper ratio (ct/cr) 

coordinates of intersection of extended leading and 
trailing edges 


local chord 


in x-direction 



S 

A 

dz/dx 

m 


wing area 

aspect ratio (t 2 /S) 
slope of airfoil surface 


slope of— source line 


yi 

Xi - £ 


) 
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mo = cot A 



t/c 

thickness ratio 


I 

source -strength factor 


u 

horizontal perturbation velocity (0 X in 

x,y- plane) 

u 

u caused by source line with reversal in 

sign of y 

U, 

M 

u caused by source line originating at point (|,t|) 

1/D 

integration operator 

. \ 


p 

pressure coefficient f -~ ^ u ) 


c d 

section wave-drag coefficient 

- 

CD 

wing wave -drag coefficient 


R.P. 

real part of complex expression 



Subscripts used with u and u indicate origin of source line. 


ANALYSIS 

General Considerations 


Within the assumptions of small disturbances, the linearized 
equation for the velocity potential is 

(1 - M 2 )/^ + 0yy + 0 ZZ = 0 (1) 

Because of the linearity of equation (l) a solution may be used to 
denote one of the •velocity components rather than the velocity potential. 
(See reference 7- ) Tlie specification of one component in this manner 
actually describes the whole flow field since the other components may 
be obtained by integrating the given component to obtain the velocity 
potential and then differentiating the results along the desired 
directions to obtain the desired components. It was shown in refer- 
ence 1 that a solution to equation (1) corresponding to a disturbance 
velocity at z = 0 of an oblique line of pressure sources in the 
plane z = 0 is, for mg >1.0, equivalent to 


I 


cos 


-1 


x - mg 2 ?- 
P(mx - y) 


u = R.P. 


( 2 ) 
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and for mp <1.0, 


where 


u 


R.P. 


I cosh -1 


x - mP^ 
P y - mxj 


I 


_ V m dz 

" ^ \t(mV - l) ** 


(3) 


(4) 


The + and - signs in equation (4) refer to equations (2) and (3), 
respectively. The origin of- the lines of sources described in equa- 
tions (2) and ( 3 ) is the origin of the coordinate system. Equa- , 
tions (2) to (4) were given in references 2 and 3 for a line of sources 
with origin located at any point ( £,q) in another form and are for the 
supersonic leading edge (m3 >1} 


qfoy) = r.p. 


1 cos' 1 x - Li ^ g (y zJL 1 


p gx - I) - (y - nj] 

and for subsonic leading edge (mp < l) 


x cosh- 1 f - J 

p|m(x - I) - (y - T))[ 


u M (x,y) = R.P. 
where for an untapered wing m = cot A and for a tapered wing 


(5) 


( 6 ) 


m 


yj 

x i - 1 


( 7 ) 


In reference 1, the following is stated: "Curved surfaces require 

a continuous distribution of sources and sinks alined with the 
generators of the surface. Each elementary source line causes an 
infinitesimal change in direction of the surface and hence the slope at 
any point may be obtained by adding up the effects of all sources ahead 
of that point." Thus, it- was shown (reference 1} that 


dz _ it f 1 V±( m g p 2 - 1) dl . g 
dx V J m d£ a 


( 8 ) 


and 


dl _ V 

d£ it 



( 9 ) 
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Pressure Distribution 
The pressure coefficient is given by 

P = - |£> (10) 

where the s umma tion in equation (10) represents the superposition of 
the finite sources of strength sufficient to form the desired angle of 
intersection of the arcs of the leading and trailing edges together with 
a continuous distribution of infinitesimal sources and sinks along the 
chord line to produce any desired curvature. 

The following equations from references 3 and 2, respectively, 
when applied to wing plan forms where each tip affects solely its own 
half of the wing can be given as follows: For a supersonic leading 

edge (moP > 1) 



^,0 + ^*0,0 " D^O “ D^jO " ^b/2mQ } -b/2 + 1^/2^, b/2 


(H) 


and for a subsonic leading edge (idqP < 1) 


y u = ^,0 + "0,0 + u c,0 + u c,0 “ ^0,0 - d“0,0 + 

D u c,0 + p u c,0 ~ u b/2mg,b/2 + r $ 1 'b/2mQ,'b/2. ^ 12 ) 

where the unbarred quantities in equations (ll) and (12) refer to con- 
tributions of the sources in the half-wing where the point under con- 
sideration is located, and the barred quantities refer to the contri- 
butions of those sources in the other half-wing, the disturbances of 
which are propagated to that point. The continuous distributions of 

sources and 6inks (represented by ^u and -j!ju) in equations (ll) and (12) 

are obtained by chordwise integration of the solution for an elementary 
line source given by the expression: 


(*,y) + U|^(x,y) = R.P. 


21 cos - ^ 


i I (* - S) g - P 2 (y - n) g 

g V m 2 (x - I) 2 - (y - t\) 2 


(13) 


The following formulas can be used to evaluate the contributions 
of the various terms of equations (11) and (12): For mgP > 1 


u o,o ( *' y) + u o,o (x,y) 


R.P. 


21 cos"^ 


1 l x 2 - P^y 2 
m^ 2 - y 2 


( 1*0 
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For MqP < 1 

+ u<),o( x >y) 



For the region ahead of the Mach line from the wing apex (x i 3y 
and diqP > 1) 



The contributions ^u + j of the chordwise distributions of source 

and sink lines are obtained by a summation of the effects of all the 
elementary source and sink lines, the disturbances of which are propagated 
to the point (x,y) . For instance, in the region between the leading edge 
and the Mach line from the wing apex (equation (l 6)) the contributions 
of the leading-edge source lines is lit, and the increment due to the 
change in source strength from an infinitesimal source line to another 
approaches it dl in the limit. In the limit, the summation becomes 

it 

the root chord. For the region behind the Mach line from the wing 
apex (x ^ j3y and ihqP > l) 




along 


an integration over the proper limits of 
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where m is a function of £ given by equation (7) • The first term 
in equation (17) is the contribution of the elementary source lines 
from the leading edge to the Mach line emanating from the wing apex, 
and the last term represents the contribution of the elementary source 
lines between the Mach line from the wing apex and the forward Mach 
line of the point (x,y). For a subsonic leading edge (ioqP < l) 


r . 


^o,o( x >y) + | i] o,o< x >y) = 


n6=x - py % cosh -i ,»-t -ifr 

di 2 ' Plm(x - |) - y 


vJ 

V 


1=0 




di 


\ 


i=o 


— — cosh 
« s 

-L. — - * 

p[m(x - |) + y_ 

\ 

' \mpj 

2 

1 - 1 


di 


(18) 


J 


where equations (l 6 ) to (l8) were given in another form in reference 1. 


In calculations for a parabolic-arc section in the stream direc- 
tion, d 2 z/d|2 is constant and may be removed from the integrands of 
equations (l6) to (18). Because of the complexity of the integrations 
(equations (l6) to (l8)) for the contribution of combined plan-form 
taper and curvature of profile to the pressure distribution, graphical 
integration was considered the most advantageous means of obtaining the 
results. 


Drag 

Since the pressure distribution (equation (10)) for the tapered 
configuration was not available in closed form, the section wave -drag 
coefficients and the wing wave -drag coefficient were calculated by 
graphical means from the equations: 
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where 


and 


§ 2 . - 

dx c\2 c J 




( 20 ) 


( 21 ) 


RESULTS AND DISCUSSION 


Pressure-distribution calculations have been made for the wing 
configuration of figure 1 having a symmetrical parabolic-arc section. 

The chordwise pressure distributions at the root- chord, tip chord, and 
four intermediate spanwise stations are presented in figure 2. The 
distribution of section wave-drag coefficient (equation (19)) is shown 
in figure 3 and was integrated graphically (equation (21)) to give the 
wing wave-drag coefficient Cp = 5.34(t/c) 2 for the wing configuration 
of figure 1 or Cp = 4.86(t/c) 2 cot A for a limited family of tapered- 
wing configurations with symmetrical parabolic -arc sections 
satisfying X = 0.531* Af5 = 4.65, and p cot A = 1.375- 

As a point of interest, the pressure distributions were calculated 
for the previously mentioned stations by an arbitrary approximation 
and are presented in figure 2j the pressure distribution at each station 
was calculated by the formulas of reference 1 for untapered wings having 
the same span and leading-edge sweep as the tapered wing and having a 
constant chord equal to the local chord of the tapered wing at the 
spanwise station under consideration. The untapered wings used in 
these computations are illustrated in figure 4. Although the results 
compare favorably in the vicinity of the root chord, the agreement 
becomes poor with distance from the root-chord. (See fig. 2.) 

The distribution of section wave-drag coefficient,- also calculated 
with the use of the geometric assumptions indicated by figure 4 and 
the formulas of reference 3* is shown in figure 3- The corresponding 
wing wave-drag coefficient is Cp = 5-29(t/c) 2 cot A * which is about 
9 percent higher than the result of this paper. 

The wing wave -drag coefficient calculated in the present paper 
(Cp = 4.86(t/c) 2 cot a) was compared with the value obtained for a 
geometrically similar wing with rhombic profile (reference 5) multiplied 
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by a profile -correct ion factor. This factor was chosen as the ratio of 
the wing wave -drag coefficients of two untapered wings (having the same 
aspect ratio and leading-edge sweep as the tapered wing) with parabolic- 
arc (reference 3) and rhombic sections (reference 6), respectively. The 
wing wave-drag coefficient Cp = 4.85(t/c) 2 cot A so obtained 
■when Af3 = 4.65 and 3 cot A = 1.375 compared excellently with the 
value calculated herein. 


CONCLUDING REMARKS 


The study of the supersonic pressure distribution and wave drag of 
sweptback wings at zero lift has been extended by means of linearized 
theory to include the combined effects of plan-form taper and curvature 
of profile. Calculations have been made for a specific wing with a 
supersonic leading edge and these results are applicable to a limited 
family of plan forms. Comparisons have been made of the pressure dis- 
tributions and the wave drag calculated by this exact linear solution 
with the results of an arbitrary simplified approximation based on 
formulas for untapered plan forms with symmetrical parabolic-arc sections. 
It appears that these approximations are adequate only in the vicinity 
of the root chord. In addition, the wing wave-drag coefficient was 
compared with a value obtained from that of a geometrically similar 
sweptback tapered plan form with a rhombic profile by using as a profile- 
correction factor the ratio of the wing wave -drag coefficients of two 
untapered wings (having the same aspect ratio and leading-edge sweep as 
the tapered wing) with parabolic-arc and rhombic sections, respectively. 
This approximation indicated excellent agreement for the cases considered. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., November 1, 194-9 
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Figure 1.- Tapered-ving configuration with supersonic leading edges; A = 3.72, A » ^2.3°, 1 = 0.531, 

and M » 1.6. 


NACA W 2009 



12 


NACA ITT 2009 




Ffercent chord 


Figure 2.- Illustrative chordwise pressure distribution on a symmetrical 
parabolic-arc section at zero lift, idq3 = 1.375* 
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